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Abstract 

We discuss the integrable boundary conditions for the one-dimensional (ID) 
Hubbard Model in the framework of the Quantum Inverse Scattering Method (QISM). 
We use the fermionic i?-matrix proposed by Olmedilla el al. to treat the twisted 
periodic boundary condition and the open boundary condition. We determine the 
most general form of the integrable twisted periodic boundary condition by consid- 
ering the symmetry matrix of the fermionic i?-matrix. To find the integrable open 
boundary condition, we shall solve the graded reflection equation, and find there 
are two diagonal solutions, which correspond to a) the boundary chemical poten- 
tial and b) the boundary magnetic field. Non-diagonal solutions are obtained using 
the symmetry matrix of the fermionic i?-matrix and the covariance property of the 
graded reflection equation. They can be interpreted as the 5*0(4) rotations of the 
diagonal solutions. 
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§1. Introduction 



The one-dimensional (ID) Hubbard model is one of the most important solvable models in 
condensed matter physics. The ground state energy of the ID Hubbard model 

N N 1 1 

ff="EE ( C L C m+ls + 4n+lsCms) + U £ (n m] - -)(n m | - -), (1.1) 
m=l a=tl m=l ~~ ~ 

with periodic boundary condition 

c N+ls = C L C N+ls = Cis (S =tl), (1.2) 

was obtained by Lieb and Wu |]TJ by means of the coordinate Bethe ansatz method. Here 
c ms and c ms are the fermionic creation and annihilation operators with spin s(=T!) at site 
m(= 1, 2, ■ • • , N), and n ms is the number density operator 

T^ms (l*^) 

The parameter U is the coupling constant describing the Coulomb interaction. The bulk prop- 
erties of the ID Hubbard model have been clarified by analyzing the associated Bethe ansatz 
equations |2J. 

As is well known, the Hamiltonian Ql. If) enjoys two SU(2) symmetries @^7|, which are the 
spin-£7£7(2) generated by 

N N i N 

m=l m=l m=l 

and the charge- SU (2) generated by 

N N -r N 

v + = Et-iritii' v~= E(- i r c mi c mT> EKt + M-i). (i.5) 

m=l m=l m=l 

For the consistency of the definition of the charge- £7 i7 (2), it is necessary to assume that the 
number of the lattice sites N is even. The spin-£7£7(2) and the charge-£7[/(2) are connected 
through the partial particle-hole transformation 



c m ^c mh - {-l) m c mi , U^-U. (1.6) 

Since a constraint 

S z + rf = integer (1.7) 



holds, the exact symmetry of the Hamiltonian Q) is 50(4) = [SU(2) x 31/(2)] /Z 2 . The 
5*0(4) symmetry may be the most fundamental property of the ID Hubbard model that char- 
acterizes the various physical features of the model. For example, it was proved by EBler et al. 
that the Bethe ansatz states of the ID Hubbard model are incomplete and have to be comple- 
mented by the 50(4) symmetry P-|10|. Efiler and Korepin fi~l"|,|i~2|| showed that the elementary 
excitations of the half-filled band constitute the multiplets of £70(4). 

To discuss the bulk properties of the model, we usually assume the periodic boundary 
condition ( |1.2p , since the bulk physical quantities do not depend on the boundary condition. 



1 



However, when we investigate the surface critical phenomena, or the transport properties, it is 
necessary to pay attention to the boundary condition. 

The twisted periodic boundary condition for the Hubbard model was discussed by Shastry 
and Sutherland [ 13L 111 . They have shown that the coordinate Bethe ansatz method is applicable 



if we generalize the periodic boundary conditions to the twisted periodic boundary condition 



,t 

■-jv+it 



e r c 



l|i c iV+lj 



-lip 



Cl|, 



;i.8) 



where the parameters <j) and ip are twist angles. These changes of the boundary condition affect 
directly the transport properties such as the stiffness constants fT^,|l4||. 

Recently the ID Hubbard model with the open boundary condition has attracted much 
attention. The corresponding Hamiltonian reads 



H, 



open 



N-l N 
^ ^ ( C lns C rn+ls + C m+ i s C ms ) + U ^ | // 



m=1 «=TI 



m=l 



1 



m T 2 ) \ n ml 2 



1 



-Pit 



PU n H 



Pni n Nl 



;i.9) 



where p\ s ,PNs (s =ti) are the boundary fields that specify the boundary conditions. In the 
case of vanishing boundary fields p\ s = pn s — (s =| j), this model was solved by Schulz JLS 



applying the coordinate Bethe ansatz method. More recently, Asakawa and Suzuki [IB] ex- 
tended the solution to the boundary chemical potential case pif = pij = p^i = Pni = P- The 
present authors 0] studied the consistency conditions for the Bethe ansatz wave functions and 
classified the possible forms of the boundary fields. It turned out that the boundary fields can 
be either of the form a) a boundary chemical potential: p m ^ = p m ^ = p m , or b) a boundary 



magnetic field: p. 



m] 



'Pull Pm: 



where m = 1 or m = N. 



Since the boundary conditions at the different ends are independent, we have four integrable 



open boundary conditions |T7 



Case A. p 1T = p n = p ± , p N] = p Nl = p N , 

Case B. p n = -p n = p u p m = -p Nl = p N , 

Case C. px T = -p H = p ls p N ^ = p Nl = p N , 

Case D. p 1T = p n = p u p N1 = -p Nl = p N . 



(1.10) 



The corresponding Bethe ansatz equations were derived Hl7|, [18[|. 

The exact integrability of the ID Hubbard model with periodic boundary condition was 



established by Shastry |T9H2T| and Olmedilla et al. [^, ^] . The Jordan- Wigner transformation 



-ml 



m—l)® 'mi 



changes the fermionic Hamiltonian ( |1.1|) into an equivalent coupled spin model 



H 



N 
m=l 



N 

E( 

m=l 



T, 



m+1 to 



T + T T 
'to 1 'm i 



m 1 m+1 / 



u N 

- V a z r z 

A / y " TO ' ml 

4 m=l 



;i.n) 



(1.12) 
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where a and r are two species of the Pauli matrices commuting among each other, and 

°"m = l^Pm ± icr m), r m = TT^m ±ir m)- (!- 13 ) 

Shastry constructed the L-operator and the i?-matrix which satisfy the Yang-Baxter relation, 
U RLL = LLR" , for the equivalent coupled spin model ( |1.12| ). The present authors proved the 
Yang-Baxter equation, U RRR = RRR\ for Shastry's LVmatrix |24H26|]. 

It is well known that the Jordan- Wigner transformation ( |1.11[ ) is not consistent with the 
periodic boundary condition for the fermion operators (|1.2|) . In fact the periodic boundary 
condition for the fermionic operators ( |1.2|) corresponds to a sector dependent twisted boundary 
condition for the coupled spin model ( |1.12|) 

f II <) 4 +1 = of, ( n r± +1 = rf. (1.14) 

\m=l / \m=l / 



Hence it is more natural to use the fermionic formulation developed by Olmedilla et al. [22|. 23]. 
when we discuss the boundary conditions for the ID Hubbard model. They proposed the 
fermionic L-operator and the fermionic .R-matrix which satisfy the graded Yang-Baxter relation. 
The transfer matrix which corresponds to the fermionic ID Hubbard model with the periodic 
boundary condition is constructed by taking the supertrace of the monodromy matrix |22],[23||. 

In this paper, as a first problem, we shall generalize their results to the case of the twisted 
periodic boundary condition. We shall find the symmetry matrix of the fermionic .R-matrix, 
which is the constant solution of the graded Yang-Baxter relation. It is closely related to 
the SO (4) symmetry of the Hamiltonian ( |1 . 1|) . The symmetry matrix can be inserted into 



the transfer matrix without breaking the integrability |27|. We shall show that the periodic 
boundary condition is twisted by the insertion of the symmetry matrix. 

A general method to prove the exact integrability of the model with open boundaries was 
developed by Sklyanin ||27|| . Through the Jordan- Wigner transformation, the Hamiltonian ( |1.9|) 
is transformed into 



N-1 N-1 jj N 

H op en — X/ (°m°ro+l + °m°m+l) + X/ ( T m T ?n+l + T m T m+l) + 7 E <J m T m 
m=l m=l m=l 

- V -f°l- P -f°»- P frl- V -fTt,. (1.15) 

Zhou |28[ applied the Sklyanin's formalism to investigate the exact integrability of the model 
(|1.15|) . He formulated the reflection equation in terms of Shastry's i?-matrix and found a 
solution which corresponds to the boundary chemical potential. Subsequently, the solution 
corresponding to the boundary magnetic field was also found [[29(1 . 

In this paper, as a second problem, we shall formulate the graded reflection equation in terms 
of the fermionic i?-matrix. It is shown explicitly that there are only two diagonal solutions for 
the graded reflection equation, which correspond to a) the boundary chemical potentials and 
b) the boundary magnetic fields. In this way, the integrability of the ID Hubbard model 
with boundary fields ( |1.9|) is proved. Moreover, making use of the covariance property of 
the reflection equation, we construct two non-diagonal solutions. It is shown that the partial 
particle-hole transformation for the fermionic i?-matrix relates these solutions. 
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§2. Fermionic .R-Matrix and Graded Yang-Baxter Relation 



Let us recall the fermionic formulation of the integrability of the ID Hubbard model f22|,|23 
The fermionic L-operator is given by 



£ m (6)=C(6)(C m i(6)®C mi (e))C(8), 



where 



and 



C{9) = cosh A/ + sinh Aa* ® a 2 



-fml(0) 



fmi(0) 



IC 



Cm J, 



(2.1) 

(2.2) 
(2.3) 



Here / is the 4x4 identity matrix, ® means the usual direct product, and 

fms(0) = sin0 - {sin9 - icos6}n ms , 
gms(9) = cos 8 — {cos 6* + isin^} n ms . 

The symbol (g> denotes the Grassmann (graded) direct product 



(2.4) 



[A® B] = (-l^M+WM^fl,,, 
P(l) = 0, P(2) = l. 
It is sometimes convenient to write the fermionic L-operator explicitly, 



(2.5) 



/ -e h fm]{9)fml{6) -fm]{0)c m l 



-i/mt(^)ii 
~ le C m1 C ml 



cU9m.l(0) 



fml (9) 



16 C m j C m ^ 



1C, 



9m\ 

' l 9m}{0)cl 



ml 



(2.6) 



rn | 



9 m] 



The parameter h should be considered as a function of the spectral parameter 9 and the 
Coulomb coupling constant U through the relation 



sinh 2h 



U 
7' 



sin 29 

The fermionic L-operator fulfills the graded Yang-Baxter relation p2 

n 12 (9 u 9 2 )[c m (9i)®c m (92)} = [c m (9 2 )®c m (9 1 )]n 12 (e 1 ,e 2 ) 

s s 

Here the parity of the Grassmann direct product ® is assigned as 

s 

P(1) = P(4) = 0, P(2)=P(3) = 1. 
In the graded Yang-Baxter relation ( |2.8|) , the constraints 

sinh 2hi sinh 2h 2 U 
sin 29 l ~ sin 29 9 . ~ 



sin 29 2 
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(2.7) 



(2.8) 



(2.9) 



(2.10) 



are assumed. The matrix elements of the fermionic R- matrix lZx 2 (9i, 2 ) are given in ref. 22. 
For later use, we introduce an equivalent fermionic i?-matrix 



n 12 {e 1 ,e 2 ) = v 12 n 12 {e 1 ,e 2 ), 

where Vu is the graded permutation 

V a ^ 5 = (-i) p ( a ) p W<UV 



(2.11) 



(2.12) 



Using the fermionic i?-matrix TZi 2 (6i,6 2 ), the graded Yang-Baxter relation ( |2.8|) can be ex- 
pressed as 

^12(^1,^2) (a*^)®/) (7®£ m (02)) = (/®£ m (0 2 )) U m {ei)®i)n 12 {e x ,e 2 ). (2.13) 

We parametrize the matrix elements of the fermionic i?-matrix as follows 

^12(01,02) = 



a+ 















































\ 





-ib+ 








e 









































-i6+ 

















e 


































-c+ 








if 








-if 








d+ 
















e 








ib~ 




















































— a~ 











































if 








c~ 








-d~ 








-if 


































i&- 

















e 
















e 

















i&~ 


































-if 








-d- 








c~ 








if 











































— a,- 




















































ib- 








e 
















d+ 








-if 








if 








-c+ 


































e 




























































e 








-ib+ 





/ 















































a+ 



where 



cos 2 ^ -6 2 ){l± tanh(/H - M^^TT 
( cos(0i - 6> 2 

sin(0i - 2 ) cos(#i - 62) jl ± tanh(/i! - h 

{COS 
1 ±tanh(/ii + h 2 ) — 
cos 

sin 2 ^ - B 2 ) (l ± tanh(/n + h 2 f^ Q l + °f " 

1 ± tanh(/ii - h 2 ) — , 
cos(0i + 9 2 ) 

, , . sin(6 l i — 60) 

sm(6'i + 6/2) 

cos(0 1 - gg) sin(gi-g 2 ) 

cosh(/ix — ft, 2 ) ' cosh(/ii + h 2 ) 




(2.14) 



(2.15) 
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The second equalities for the Boltzmann weights b ± and d ± are due to the constraints (|2.10| ), 
or equivalently the relation 

tanhQi - h 2 ) _ tan(0i - 2 ) 
tanh(/ii + h 2 ) ~ tan(0i + 9 2 ) ' 

We note some useful relations among the Boltzmann weights [p^j , 

d± = a ± + c ± , d + d~ = e 2 + f, 

+ b + b~, f = b + b- + C+C". (2.17) 



2 + - 

e = a a 



The monodromy matrix is defined as the ordered product of the fermionic L-operators 

N 

T(0) = n C m {9) = C N {9) ■■■C l {9). (2.18) 

m=l 

From the (local) graded Yang-Baxter relation ( |2.8| ), we have the global relation for the mon- 
odromy matrix 

£i 2 (0i,02)[T(0i)®r(0 2 )] = [T(9 2 )^T(9 1 )}n 12 (9 u 9 2 ), (2.19) 

s s 

or equivalently 

n 12 (9i, 9 2 ) T (9t) T (9 2 ) =T (9 2 ) T (9 1 )TZ 12 (9 1 , 9 2 ), (2.20) 

where 

T (9) = T(9i) ® J, T (0) = I ®T{9 2 ). (2.21) 

s s 

In other words, the monodromy matrix ( [2.18 ) is a representation of the associative algebra T 
defined by ( pA9|) or ( ggO|) . 

By taking the supertrace of ( [2.19D or (|2.20|) , we find that the transfer matrix 

t{9) = strT(0) = tr {(a 2 ® a") T(0)} (2.22) 

constitutes a commuting family 

[t(0i),t(0 2 )] =0. (2.23) 

Olmedilla et al. showed that the Hamiltonian ( |1 . 1|) under the periodic boundary condition 
( |1 . 2|) can be obtained by the series expansion of the transfer matrix t(9) around 9 = 0. 
Now we list some properties enjoyed by the fermionic i?-matrix TZ\ 2 (9i, 9 2 ). 

(1) Regularity (Initial condition): 

n 12 {9 ,9 )=V 12 . (2.24) 

(2) Graded Yang-Baxter equation: 

n 12 (9r, 9 2 )TZ 13 (0i, 9 3 )TZ 23 (9 2 , 9 3 ) = TZ 23 (9 2 , 3 )ft 13 (0i, 3 )^i 2 (0i, 2 ), (2.25) 

where we assume the constraints 

sinh2/ii sinh2/i2 sinh2/i3 U 



sin 20x sin 20 2 sin20 3 4 
6 



(2.26) 



(3) Unitarity: 

where 

and 



p{9t, 9 2 ) = cos 2 (#i - 9 2 ){ cos 2 (9 1 - 9 2 ) - tanh 2 (/n - h 2 ) cos 2 (9 1 + 9 2 )}. 
(4) Crossing unitarity: 

nfi{-9 u 9 2 )^{9 2 ,-9 1 ) = nfi(-9 u e 2 )-R%{p 2 , -d 1 ) = p(e 1 ,e 2 ) i, 

where 

p(0 u 9 2 ) = sin 2 (^ + 9 2 ){ sin 2 (^ + 9 2 ) - tanh 2 (/n - h 2 ) sin 2 (^ - 9 2 )}. 
Here the supertransposition acts on the 4x4 matrix with the parity Q2.9p as follows 

/ Dn C\\ C\ 2 D\ 2 \ 

B\\ A U A\ 2 B 12 

B 2 \ A 2 \ A 22 B 22 
V D 2 i G 2 \ C 22 D 22 J 

In ( p.30| ), stj means the supertransposition with respect to the j-th space. 
We further note some properties of the .R-matrix, 





-B u 


—B 2 \ 


D 21 \ 


C n 


An 


A 21 


C 2 \ 


C\ 2 


A 12 


A 22 


c 22 


V D 12 


—Bi 2 


—B 22 


D 22 J 



K 12 (9 1 + ^9 2 + ^] 

- 9 -- 

2' 2 2' 



Tl 12 {9n9 2 ) = 

= n l2 {9, 

H 12 {9n9 2 ) st ^ = fti 2 (#iA) stl ' st2 . 
Here stj stands for the inverse of the supertransposition stj, 

(x st ) st = (x^) st = X. 



(2.27) 
(2.28) 
(2.29) 

(2.30) 
(2.31) 



(2.32) 



(2.33) 
(2.34) 

(2.35) 



We remark that TZi 2 (9i, 9 2 ) stl,st2 can be obtained from the .R-matrix TZi 2 (9i, 9 2 ) by exchanging 
the Boltzmann weights as if <-> —if. 

The graded tensor product is assumed in the above relations. The graded Yang-Baxter 
equation, for instance, is expressed as follows [SO] 



Ka-MOu 2 )n alc , jd {9 u 9 3 )K blc ,., kl (9 2 , 3 )(-)^ 6 ')(^)+iV)) 
= K bc * ld {9 2 , 9z)K^ an {9 u 9*)Ka>v>Ah, 9 2 )(-) p ^ p ^ +p ^\ 



(2.36) 



The authors gave a proof of the Yang-Baxter equation for Shastry's .R-matrix p4|-Pq1 . We can 
show that the graded Yang-Baxter equation for the fermionic .R-matrix is equivalent to the 
Yang-Baxter equation for Shastry's .R-matrix. 
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§3. Symmetry of the Fermionic .R-Matrix and Twisted Periodic 
Boundary Condition 



In this section we shall discuss two important symmetries of the fermionic i?-matrix. The 
symmetry of the fermionic i?-matrix is defined by the constant matrix M = (My) that satisfies 



n l 2(e 1 ,e 2 ),M®M 



0. 



or equivalent ly, 



n l2 {9 1 ,9 2 ){M®I){I®M) = (I®M)(M®I)K 12 (6 1 ,6 2 ) 



(3.1) 



(3.2) 



Here we assume that the matrix elements My are commuting numbers. Solving the defining 
relation ( |3.1| ) or ( |3.2| ), we have found that the symmetry matrix M for the fermionic i?-matrix 
takes the following form 

f Mu M 14 \ 
M 22 M 23 
M 32 M 33 
V M 41 M 44 J 



M 



(3.3) 



with the condition 

M11M44 - M41M14 = M 22 M 33 - M 23 M 32 . 
We denote the submatrices of M as 



(3.4) 



M 



11 



charge 



Mn M; 
M41 M44 



M 



j ^"spin 



M 22 M 23 
M 32 M 33 



Then the condition 



can be written 



detM, 



charge 



detM 



spin 



A(M). 



Since an overall constant is not relevant for the integrability, we put 

A(M) = 1. 



(3.5) 



(3.6) 



(3.7) 



In other words, the submatrices M c h ar ge, M spin belong to SX(2,C). The symmetry matrix M 
reflects the SO (4) symmetry of the Hamiltonian. It is to be remarked that the symmetry 
matrix M of Shastry's i?-matrix Ri 2 (9\, 8 2 ), which is defined by 



R 12 (9 1 ,9 2 ),M®M 



0. 



(3.8) 



is not of the form ( p.3|) and (|3.4j ) plfl . Shastry's i?-matrix Rx 2 {9\, 9 2 ) is related to the fermionic 
i?-matrix 7^-12(^1, ^2) through the formula p2| 



Ri2(9i,9 2 ) = W 1 ~ 2 1 n 12 (9 1 ,9 2 )W 1 



12: 



(3.9) 
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where W± 2 is a diagonal 16 x 16 matrix 

W 12 = diag(l, 1, -i, -i, -i, -i, 1, 1, -1, -1, i, i, i, i, -1, -1). (3.10) 

Since ( |3.9| ) is not a gauge transformation, the symmetry of Shastry's i?-matrix may be differ- 
ent from that of the fermionic .R-matrix. The fermionic i?-matrix is more suitable when we 
investigate the symmetry of the fermionic Hamiltonian in the framework of the QISM . 

Using the symmetry matrix, one may twist the periodic boundary condition. The symmetry 
matrix M can be inserted into the transfer matrix without breaking the integrability \27\ 

t(6; M) = str {MT(6)}. (3.11) 

The generalized transfer matrix t(9; M) still constitutes a commuting family 

[t(6 1 ;M),t(e 2 ;M)]=0. (3.12) 

The series expansion of t(6; M) around 8 = gives rise to a Hamiltonian H(M) 

t(6; M) = str {MT(0)} (l + 9H(M) + • • •) (3.13) 

where 

N-i N 1 1 

H{M) = - ( c Lc m +is + c ] m+ls c ms ) ( n mt - ^){n m \ - -) 

m=l s=U m=l A A 



- [M n M 22 c\^ + iM 11 M 32 c t u - M 41 M 32 c lt - iM 41 M 22 c u ) c m 
-c ] N} (M 33 M 44 c 1T + iM M M 23 c n - M 14 M 23 cl T - \M U M&^ 

- (MnM^c^ - iM u M 23 4 T - M 4 iM 23 c u + iM 41 M 33 c u ) c m 

-4; (M 22 M 44 c u - iM 44 M 32 c 1T - M 14 M 32 c f u + iM 14 M 22 c 1T ) . (3.14) 

It is easy to see that the choice M = I corresponds to the periodic boundary condition (|1.2p 



Further, we assume the submatrices M c h ar ge and M spin belong to SU{2). Namely, 



M 44 


= M* u , 


M 41 = 




\M 1X 


2 + 


|M 14 | 


M 33 


= M* 2 , 


M 32 = 


-m; 3 , 


\M 22 


2 + 


|M 23 | 



2 



2 



1, (3.15) 
1. (3.16) 



Here the symbol * means the complex conjugation. In this case, we find that the periodic 
boundary condition is twisted as follows 

4 +1T = MnM 22C ; T - iMnM^Si - M* 4 M 2 * 3 c 1T + iM* 4 M 22 c u , 

Cjv+ i T = \r A \ I ■:■><■.■ + iM^Mnexx - M 14 M 23 c\i - iMuM*^, 

4 +u = M 11 M; 2 c\ l - iM n M 2 ^ + M* 4 M 23 c u - iM* 4 M* 2 c lt , 

c N+n = M 22 M* lCll + iM* 1 M* 3 ciT + M 14 M*3C f u + iM 14 M 22 4 r (3.17) 
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The relations Q3.17 ) can be expressed in a matrix form 



cjv +u iCAr+lt 

t 



M 



1C 



JV+lt c N+ll 



spin 



C 

ic 



h 1c it M 

t ^ I charge- 



(3.18) 



This shows that the periodic boundary condition can be rotated by the group SU{2) x SU{2). 
A quite similar SO (4) rotation was observed by Affleck 0. 

Since the choices M spin = — l,M charge = 1 and M spin = l,M charge = — 1 induce the same 
transformation fl3.18| ), the exact group symmetry is SO(4) = [SU(2) x SU(2)]/Z 2 . 

Now let us consider the diagonal M parametrized as follows, 



M = S a * <g> e 1 


2 a 








/ ; <t>+i> 

e 2 








\ 







e 2 
















e 2 







I o 








■.<t>+i> 

e' 1 — J 



(3.19) 



From ( |3.17p , the corresponding twisted periodic boundary condition is 



J 



e^c 



e^c^, cjv+ij. 



which is identical to (ITT 



Especially, if we take 



(3.20) 



it and ■?/> = —it, (|3.19|) becomes 



M 



{ 1 \ 
0-100 
0-10 

V o o 01/ 



(3.21) 



In this case, the supertrace in the transfer matrix ( |3.11| ) becomes the usual trace 

str{(cr 2 <g> a z )T{6)} = tiT(0). 



(3.22) 



Thus, as remarked in ref. 23, we can take the trace of the monodromy matrix instead of the 
supertrace. However, the boundary condition in this case should be anti-periodic 



-N+ll 



— c 



C 1T> C N+11 

t 

14.) Cn+ii 



-Cit, 
-c U . 



(3.23) 



Next we consider a discrete symmetry related to the partial particle-hole transformation 
TT6|) . We solve the following equation 



n 12 (6 1 ,6 2 ;U)(N®I)(I®N) = (I ®N)(N ®I)K 12 (6 1 ,8 2 ;-U). 



(3.24) 



with a constant matrix iV = (Nij). Here we explicitly write the [/-dependmence of the fermionic 
i?-matrix. Note that the coupling constant of the fermionic i?-matrix in the RHS is —U, or 
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equivalently hi — > — h±,h 2 — > — ^2- By solving the denning realtion (|3.24j ), we find that the 
constant matrix N has a form, 



N 



( 





iVl 2 


iVl3 





\ 


iV 2 i 








iV 24 






N 31 












\ 





iV 42 


iV 43 








with a condition 

We set 



A(N) = N 12 N 43 - N 13 N 42 = N 2 iN u - N 31 N 2i . 

1. 



A(N) 

as before. The matrix N can be written as 



(3.25) 

(3.26) 
(3.27) 



N 



a 



( N 21 








N 24 


\ 





N 12 


N 13 








N 42 


iV 43 







V N 31 








A^ 3 4 


/ 



(3.28) 



Here and hereafter, 1 means the 2x2 identity matrix. Thus N is a composition of the symmetry 
matrix ( |3.3| ) and 

^0100^ 

Q = Q- 1 = l 



a 



(3.29) 



10 
1 
V 1 J 

In terms of this matrix, the fermionic i?-matrix with the coupling constant U is transformed 
into the fermionic i?-matrix with the coupling constant — U as follows 



Q-^KniPu 9 2 ; U) QQ= n 12 {9 1 , 9 2 ; -U). 



(3.30) 



The matrix Q is related to the partial particle-hole transformation ( |1.6| ) of the Hamiltonian. 
There are other choices for the matrix Q. For example, 



Q = a z 



a" 



( -i 0\ 

i 

i 

V -i J 



(3.31) 



also satisfies (|3.30|) , which may be more appropriate as the partial particle-hole transformation 
( Of ) P2"] . However, for simplicity, we choose the matrix ( |3.29| ) as the partial particle-hole 
transformation of the fermionic i?-matrix. 
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>4. Graded Reflection Equations for the Fermionic i^-Matrix 



In this section we investigate the integrability of the ID Hubbard model with open boundary 
condition in terms of the fermionic i?-matrix 7^12(^1,^2)- We introduce an associative algebra 
71 defined by the fermionic R- matrix |27|, |33|, |34|| , 

1 



h)n 12 {e x , -6 2 )<r-(6i)n 21 (-6 2 , -e x ) 



2 

T- 



where 



1 

71 



2 

71 



(4.1) 
(4.2) 



The relation ( |4.1| ) is called the graded reflection equation (graded RE). 

The following theorem is fundamental for the application of the associative algebra 71. 
Theorem 



34 



Let 71(0) be some representation of the associative algebra 71 ( |4.1| ) and T(0) of the associative 
algebra T ( gj^ ). Then 71(0) defined by 



71(0) =T(0)71(0)T- 1 ( 



(4.3) 



is also a representation of 71 provided that the matrix elements of 71(0) and T(0) commute. 
This property is sometimes called the covariance property of the (graded) reflection equation 



34|. For the application to the ID Hubbard model, we choose T{9) and 71(0) as 



T(0) = C N (6)C N - 1 (e)---C 1 (6), 
71(0) = K.(0). 



(4.4) 



Here K_{6) is a constant supermatrix satisfying the graded RE (4.1), 



n 12 {e u e 2 )K-{e l )n 2 x{e2, -b x )k-(9 2 ) 
= K-(02)n 12 (0u -0 2 )iL(0i)7M-0 2 , -0i 



(4.5) 



The J\~-matrix K_ (0) specifies the integrable boundary condition. In this paper we assume that 
the matrix elements of (0) are commuting numbers and only the even elements with respect 
to the parity (|2.9 ) are non zero [3*5[, 



K-{6) 



( K u (6) K 14 (6) \ 

K 22 (6) K 23 (6) 

K 32 (9) K 33 (9) 

V Ku{6) Ji4 4 (0) / 



(4.6) 



To construct a transfer matrix, we introduce another J\~-matrix K + (6) satisfying the conjugated 
graded RE 



^^(02,00 (0 1 )^f 2 1 ' st2 (-0 1 ,0 2 ) Kf 2 (0 2 ) 
2 _ 1 



~-Kf 2 (0 2 )^ t 1 1 ' st2 (-0 2 ,0i) K* 1 (0i)ft 



Sti,st2 1 
12 I 



h,-o 2 



(4.7) 
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). Then the supertransposition to 



We assume that the K- matrix K + (6) is also of the form 
the i^-matrix K + (9) reduce to the usual transposition. 
Now we define the transfer matrix by 

t(6) = sir {K + (6)T46)} 

= ste{K+{6)T{6)K_{6)T- x {-0 

Using the unitarity ( |2.27| ), the crossing unitarity ( |2.30|) and the graded REs ( |4.1| ) 
show that the transfer matrix r(9) constitutes a commutative family p7|, i.e., 



[t(9 1 ),t(6 2 )}=0, 



(4.8) 

|), we can 

(4.9) 



which shows the existence of the mutually commuting conserved currents including the Hamil- 
tonian. 



We note 



where 



T 



-i. 



zr 1 

N 



(4.10) 



C~ L ( 



( - 



1 



COS 4 # 



-\c 



le C m1 C ml 



— e C to| C «4 



1 



cos 4 # 
Here we define 



c{e)[c m ^e)®c ml {B))c{e). 



QtnM 



C rnlJ m l\ 



— e 



9m]{®) C ml 



— 1C 



m 1 



iG C-mf ^mj. 
-p h n* 



(4.11) 



(4.12) 



§5. Integrability of the ID Hubbard Model with Boundary Fields 



In this section we solve the graded REs (|4.5|) and (|4.7|) for K_{6) and K + {6). First, we note a 
relation between the i?-matrix 



Rvi{Qi, #2) = R* 2 i{—9 2 -, —0i)- 
Then the graded RE (|4.5[ ) for K_(9) is cast into a form 



(5.1) 



K-{e 2 )n 



12^1; 



-e 2 )K-(e 1 )iv 12 (e 1 ,e 2 ), 



(5.2) 
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with the graded tensor products, 



K-( 



K- 



(5.3) 



Assuming K_(9) in a diagonal form 

K-{9) = I 
V 



\ 

x 2 {9) 

x 3 (9) 

x 4 (0) J 



(5.4) 



and substituting it into the graded RE ( |5.2| ), we have 10 non-trivial functional equations for 
Xi(9), (i = 1, • • -4). By solving these functional equations (see Appendix), we obtain the fol- 
lowing two sets of solutions for the diagonal K_{9). 
a) K_(9) = K {a) (9;p-): 



xi(6) 
x 2 (B) 
x 4 (9) 

b) K46) = K {b) (9;p_): 

x 2 (e) 

x 3 (6) 



(l - p^e' 2h tan 9) (l - p.e 2h tan 9 



-2h 



x 3 (6) 

(l + p_e 2/l tanfl) (l + p_e~ 2h tan9 



1 + p_e 2ft tan 9 ) ( 1 - p^e 2h tan 9 



(5.5) 



xa{9) 



■2h 



(l + p_e~ 2/l tan#) (l -p_ 



-2/i 



tan # 



1 - p_e 2h tan I 



1 + j9_e" 2?i tan ) ( 1 + p.e™ tan ) . 



1 -p_e 2/l tan6> 



2h . 



(5.6) 



Here p- is a constant parameter corresponding to the boundary field. Recall that the parameter 
h is regarded as a function of the spectral parameter 9 through the constraint 



sinh 2h _ U 
sin 29 ~ 7' 

When we put p_ = 0, the solutions K_ (9;pJ) and K_ (9;p-) are proportional to C 2 (9), 



(5.7) 



K w (9-p_ = 0) 



K y _\9-p. = 0) 



c-' i £ 2 (^), 

e' 1 /: 2 ^), 



(5.8) 
(5.9) 



which supply the missing interaction term at the boundary site m = 1. To see the contributions 
of the boundary fields clearly, we redefine 



K 

k 



M (9;p_) = e h C- 1 (9)K^(9;p_)C- 1 (9) 1 
{b) '9 ] p_) = e- h C- 1 (9)K {p) (9;p_)C- 1 (9). 



(5.10) 
(5.11) 



Expanding K ( l\9;p_) around 9 = 0, we have 



K 



l-p-(a z ®l+l®ff 2 )e + 



(5.12) 
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Now the linear term with respect to the spectral parameter produces the boundary chemical 
potential at the site m = 1 

-p_(n 1T +n u -l). (5.13) 

Similarly, we have 



K^(9;p^) = 1® l-p_((7 z ® l-l®a z )6 + ---. 
The linear term in (5.14) gives the boundary magnetic field at the site m = 1 

-p_(ni T -n n ). 



(5.14) 



(5.15) 



Thus we have shown that (9;p_) corresponds to the boundary chemical potential, while 
K ( l\9;p-) corresponds to the boundary magnetic field. 

A comment is in order on the derivations of (|5.13 ) and ([5.15 ). Rigorously speaking, ( 5.12|) 
and Q5.14D respectively give twice of (|5.13|) and G5.15|) . However, the bulk Hamiltonian appears 
twice in the expansion of the transfer matrix. Therefore the Hamiltonian contains the boundary 
fields as given in ( |5.13|) and (|5.15| ). 

We also remark that the diagonal solutions ( |5.5| ) and ( |5.6| ) coincide with those found for 
the coupled spin model ( |1.15|) |28|, p9f . Actually, if we formulate the reflection equation in 
terms of Shastry's i?-matrix [E8J, we obtain the same functional equations as ( |A.5| )-( TA.14| ) (see 
Appendix) for the diagonal K_{9). 

Solving the graded RE ( [4.7] ) for the diagonal K + (9) in a similar way, we obtain two solutions 
for the .fT-matrix K + (9) as follow, 
a) K + (d) = K[ a) [ 



xi{6) 
x 2 {9) 
x 4 (9) 



(p + + e 2h tan 0) (p+ + e~ 2h tan 9 



x 3 (9) 



2h 



(P+ + 



-21, 



p + - e 2h tan 9 ) [p+ - e~' zn tan 9 



tan^ 



-2h . 



p+- e 2h tan 



(5.16) 



b) K + (9) = K i t\9 ]P+ 



xi(9) 
x 2 {9) 
x 3 {9) 



-2h 



Xi{9) = e 
p + + e 2h tan 
p + - e 2h tan 9 



p + + e 2h tan 9) (p + 
P+ + e 

P+ - c 



3 2h tan9 



9) (p+ + e- 2h tmifl 



2h tan9 



(5.17) 



In ( |5.16| ) and ( |5.17| ), the parameter p + represents the strength of the boundary field at the 



site m = N. The solution K+\9;p + ) corresponds to the boundary chemical potential and 
(9;p + ) corresponds to the boundary magnetic field. 



As in the case of ( 5.10 ) and Q5.ll ), we may redefine 

K?(9;p + ) = e- h C(9)Ki a) (9;p + )C(9), 
K?(9;p + ) = e h C(9)K^(9;p + )C(9). 



(5.18) 
(5.19) 
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Then around 9 = 0, we have 



p 2 + (1 ® 1) + p+ ((J 2 ® 1 + 1 <g> (T 2 ) 



pi(l®l)+p+ (a 2 



1 - 1 <g> a 2 ) - (a 2 
From ( |5.2(J| ) and (3.21), we find the following properties 

stri^V = 0;p + ) = stiK { +\e = 0;p + ) = 0, 



a z ®a z )9 2 + -- 
a z ) 9 2 + -- 



d fAa), 



strl—K^(9;p 



8=0 



Btr (^ 6) (0;p + ) 



9=0 



0. 



(5.20) 
(5.21) 



(5.22) 
(5.23) 



Using the solutions (|5.5| ), ( |5.6|) , ( |5.16|) and (|5.17|) , we can construct the four integrable 
transfer matrices r(0) = str {K+{6)T{6)K-(0)T-\-6)}. 



Case A. K-(B) = K 

Case B. K_(0) = K 

Case C. K^(9)=K 

Case D. = K 



(*>)/ 
(&)/ 

(a), 



;pij 
;pi) 
;p0 
;pi) 



K + (6) = K { :\e;p N ). 



K + {6) = K ( i\6 ] p N ). 



K + {e) = K^\e ] p N ). 



K + (9) = Kf\6;p N ). 



(5.24) 



These results are consistent with the coordinate Bethe ansatz for the cases ( |1.10|) . The Hamilto- 
nian (|1.9|) is obtained by expanding the transfer matrix with respect to the spectral parameter 
9 

(5.25) 



T {9) = G x 9 + C 2 9 2 + C 3 [H opcn + const.) 9 6 + ■ ■ ■ . 

Here Cj (i = 1, 2, 3, • ■ ■) are some scalar functions. We omit the detailed derivation of (|5.25|) in 
this paper, which can be carried out using the method in ref. 22. Due to the properties ( 5.22p 
and (|5.23| ), the Hamiltonian H open appears at the third order of the expansion. 
The commutativity of the transfer matrix 



[r(9 1 ),r(9 2 )]=0 



(5.26) 



ensures the existence of an infinite number of the conserved currents in involution. The exact 
integrability of the ID Hubbard model with the boundary fields ( |1.9| ) is established in this way. 
The fT-matrices (|5.16[) and ( |5.17|) are slightly different from those for Shastry's .R-matrix 



2"8", 2"5|, which we denote K + (9). The former can be obtained from the latter by multiplying 
the diagonal matrix o z ® a z , 

K + (9) = (a z ®a z )K + (9). (5.27) 

Due to the definition of the supertrace Q2.22|) , both formulations provide the same transfer 
matrix 

str {K + (9)T(9)K49)T-\-9)} = tr {k + (9)T(9)K49)T~ 1 (-9)} . (5.28) 
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>6. Non-Diagonal Solutions of the Graded Reflection Equation 



We can construct non-diagonal solutions of the graded RE ( |4.5| ) by use of the covariance 
property (|4.3|) and the symmetry matrix M. For example, we have a solution 



K 



(«), 



where 



K n {0) 
K 22 (6) 



;P-,M charge ) = MK^{B-pJ)M- 1 

( K X1 {9) K u {9) \ 

K 22 {6) 

K 33 (6) 

V #41 (0) K u (6) J 



1 - 2p_(M n M 44 + M 14 M 4X ) cosh2/itan6> + p 2 _ tan 2 9, 

K 33 {9) = e' 2h (l - p 2 _e Ah tan 2 0) , 

1 + 2j9_(MnM 44 + M14M41) cosh 2h tan 9 + p 2 _ tan 2 9, 

4p_M u M 14 cosh 2h tan 0, 

-4p_M 4 iM 44 cosh 2h tan 0. 



(6.1) 



(6.2) 



Note that K ( l\9; p_, M c h ar ge) does not depend on the submatrix M spin . 
Since the relation 



M n M 4 4 + M 14 M 41 = { (M n M 44 - M 14 M 41 ) 2 + 4M n M 14 M 41 M 44 } 2 

= (1 +4MnMi 4 M 41 M 44 )5 (6.3) 

holds, the solution ( p.2\) depends on three arbitrary parameters, 

p_, a = M U M 14 , /? = -M 41 M 44 . (6.4) 

The corresponding boundary term is 

-p_{(M n M 44 + M 14 M 41 ) (n 1T + n u - 1) - 2iM 11 M 14 c t 1T 4 | - 2iM 41 M 44 c u Ci T } 

= -p-{(l - Aa(3)i (m T + n u - 1) - 2io4 T c T u + 2i/?c u ci T }. (6.5) 

Note the existence of the term c\+c\, (ci^cif) which creates (annihilates) a double occupied state 
at the boundary site. If we assume M c h ar ge £ SU{2) (cf. ( |3.15|) ), we have 



P = a\ 

and the boundary term ( |6.5| ) becomes hermitian (p_ is assumed to be real). 
Similarly, we have 



(6.6) 



K 



(&), 



spm / 



M#i^(6';p_)M" 1 

/ #11 (0) \ 

# 22 (9) K 23 {9) 

# 32 (0) K 33 {9) 

V K u (9) J 



(6.7) 
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where 



K 22 (9) 
K 33 (9) 
K 23 {9) 
K 32 {9) 



K M {6) 

1 - 2p_(M 22 M 33 + M 23 M 32 ) cosh 2/i tan + p 2 . tan 2 9, 
1 + 2p_ (M 22 M 33 + M 23 M 32 ) cosh 2/i tan + p 2 . tan 2 9, 
Ap^M 22 M 23 cosh 2/i tan 9, 
-4p„M 32 M 33 cosh 2/i tan 9. 



(6.8) 



Now i£l (0;p_, M spin ) does not depend on the submatrix M charge . Because of the relation 

M 22 M 33 + M 23 M 32 = (1 + 4M 22 M 23 M 32 M 33 )^ , (6.9) 
the solution ( |6.8| ) depends on three arbitrary parameters, 

p_, « = M 22 M 23 , (3 = -M 32 M 33 . (6.10) 
The corresponding boundary term is 

-p_{(M 22 M 33 + M 23 M 32 ) (m T - n u ) + 2iM 22 M 23 c t lt c u + 2iM 32 M 33 cl iClT } 

= -p_{(l - 4a/3)5 (ni T - n u ) + 2iac f 1T c u - 2i/5c f u d T }. (6.11) 

This time we see the terms cLci| and <^li c iT which hip the spin of an electron at the boundary 
site. If we assume M spin e SU(2) (cf. ( |3.16|) ), then (3 = a* and the boundary term ( |6.11|) 
becomes hermitian. 

The obtained results can be explained by considering the SO (A) rotation at the boundary 
site m = 1, 



M" 1 

spin 



C U 1C 1T 
IT C U 



M, 



ic 



charge > 



(6.12) 



where we assume M charge , M spin G SU{2). The transformation ( |6.12| ) changes the boundary 
chemical potential into 



n iT + n u 



M 



|M- 



(m T + n u - 1) - 2iM n M 14 c T 1T c T u + 2iM* 1 M* 4 c u c 1T , (6.13) 
while the boundary magnetic field is changed into 

n X | - n u = (|M 22 | 2 - |M 23 | 2 ) (n lt - n n ) + 2iM 22 M 23 c t 1T c u - 2iM* 2 M* 3 c t u ci T . (6.14) 



We see that fl6.13| ) and (6.14) coincide Q£>.5|) and ( |6.11| ) respectively. One may check that the 
Coulomb interaction term is invariant under the transformation ( |6.12| ), 



nil 



(6.15) 



as it should be. 

Thus we have obtained two non-diagonal solutions of the graded RE (|4.5|) . We have con- 
firmed that under the condition ( |4.6| ), the solutions ( |6.2[ ) and Q6.8| ) provide the most general 
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solutions of the graded RE (|4.5|) . It is also possible to construct the non-diagonal solutions of 
the conjugated graded RE ( |4.7|) in a similar way. 

Next we clarify the relationship between the two solutions K < f'\9; p_, M charge ) (|6.1|) and 
K ( l\6;p- } M spin ) ( |6.7p . Applying the partial particle-hole transformation of the fermionic R- 
matrix to the graded RE ( |4.5|) , we find that for a given solution K_(8; h(6)), 

K_{6, h{8)) = QK_(e, -h(0))Q-\ (6.16) 

is also a solution of the graded RE ( [4.51 ). Here we have explicitly written the parameter 
h{9). The transformation h(6) — > —h(9) corresponds to U — > —U. Then the solution ( |6.7| ) is 
connected with the solution (|6.1|) through the following formula, 

K^\e 1 h(e) ]P .,M spin ) = QK^\e,-h(e) ] p. J M chargc )Q- 1 . (6.17) 

Here the matrix elements of M c h arg e are exchanged with those of M sp - m . In fact one can see that 
the partial particle- hole transformation ( |1.6| ) changes the boundary term ( |6.5| ) into ( |6.11 ) with 
the exchange M charge <-> M spin . 



§7. Conclusion 

In this paper we have studied the integrable boundary conditions for the ID Hubbard model 
from the point of view of the Quantum Inverse Scattering Method. We have treated both the 
twisted periodic boundary condition and the open boundary condition. The most important 
object in the investigation is the fermionic i?-matrix for the ID Hubbard model found by 
Olmedilla et al. ||22|| . It has an interesting symmetry matrix, which reflects the 50(4) symmetry 
of the Hamiltonian. Using the symmetry matrix, we have found the general twisted periodic 
boundary condition. In a sense, the periodic boundary condition can be twisted by applying 
the 50(4) rotation to the boundary operators. We have also found the discrete symmetry 
of the fermionic i?-matrix, which corresponds to the partial particle-hole transformation for 
the Hamiltonian. Recently the 50(4) symmetry of the transfer matrix was investigated by 
Gohmann and Murakami |32|| . The symmetry of the fermionic i?-matrix is directly connected 



with the 50(4) symmetry of the transfer matrix |31 



For the integrable open boundary conditions, we have formulated the graded reflection 
equation in terms of the fermionic i?-matrix. By solving directly the functional equations, 
we have obtained the diagonal i^-matrices. There are two types of the .ff -matrices, which 
correspond to a) the boundary chemical potential and b) the boundary magnetic field. Moreover 
we have obtained non-diagonal i^-matrices using the covariance property of the graded reflection 
equation. We can rotate the diagonal boundary fields by means of the 50(4) symmetry. The 
two types of solutions are related through a partial particle-hole transformation for the fermionic 
i?-matrix. 

Recently there was reported some important progress in the evaluation of the eigenvalues 
of the transfer matrices of the ID Hubbard model |3l|[36|]. It is an interesting problem to 
generalize these results to the model with the twisted periodic boundary conditions and the 
open boundary condition. 
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Appendix: Diagonal Solutions of the Graded Reflection Equation 

We solve the graded RE ( |5.2| ) for the diagonal K^{6) 



n 12 {e u e 2 )K-{9i)K* 12 {e u -e 2 )K-{e2) 



K-{d 2 )K 



Let us introduce the Boltzmann weights of 1Z 



m<7i, 



-0 2 ) 



t ~ a+ 
o 















V o 



where 



-02) 

o 



e 






















e 














-5+ 



if 



-if 



d+ 










ib~ 


















if 



c~ 



-dr 
o 



-if 














i&" 





e 







e 






i&" 














-if 



-J- 




c — 



if 


















ib~ 


e 








d+ 




-if 



if 



-c+ 













e 






-iS+ 






cos 2 (^ + B 2 ) {l ± tanh(/ il + h 2 )^—^ 
{ cos(6>i + 0; 

wx(9i + 6 2 ) cos(9 1 + B 2 ) ( 1 ± tanh(/i! + h f 1 ^ 1 ~^ 

[ sm(t7i + v 2 

■ ( n , n\ / /i i /i \ fi i j w, , ^Os(#i -0 2 

sin(&i + U 2 ) cos 



> 1 + g 2 ){l±tanh(/. 1 -M COS ^ 1 

{ cos(6>i + 9 2 

dn 2( 0l + 2 ) | x ± tanh ( /il _ h2 f in ^ ~ ^ 1 



eP = 1 ± tanh(/ii + /i 2 



cos(^i + 9 2 



COS(#i - (7 2 J 



1+f v,^ , , sin(^ + g 2 ) 

1 ± tanh(/ii - /i 2 )^ — 77^ jtt, 

sin(fi - 9 2 ) 

cos(^i + fl 2 ) ^ _ sin(0i + 9 2 

cosh(/i 1 + h 2 



h + gg) ~ f = 

h, -u 7i„V cosh(/i! — /i 2 ) 














e 



-iS+ 





(A.l) 



\ 

o 



(A.2) 



(A.3) 
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We note useful relations among the Boltzmann weights of i?-matrix, 
,± 



e 

f 



./ 

e 

/ 
&± 

e 

&± 

7 



e ±(hi-/w) CQS ^ cog q 2 + e T(fti-/w) sin ^ sin q 2j 

e ±(h 1+ h 2 ) gin ^ cog ^ _ e T(h 1+ h 2 ) cog Qi gin ^ 
e ±(fti+h2) cos 1 cos 2 _ e T(^i+^2) gin Q x gin ^ 

e ±(fci-ha) ^n X CO s 2 - e^^ 1 -^) cos X sin 2 = sin(0i - 2 ) sin(0i + 2 



'7' 



e ±(hl+/i2) sin 6 1 cos 2 + e T{hl+h2) cos X sin 2 = sin(0i - 2 ) sin(0i + 9 2 ) 



e 

/' 

e ±(^i-ft2) cos ^ cos0 2 -e T(ftl -^sin0a sin0 2 = cos(0i - 2 ) cos(0i + 2 ) — . 

e 



(A.4) 



The graded RE ( |A.1| ) for the diagonal K^(6) is equivalent to the following set of 10 equations: 

(A.5) 



xi(0 2 ) (^x 1 {9 1 ) + yS 2 (0i)J = x 2 {6 2 ) ( b jx 1 (9 1 ) + yx 2 (0 x ) 



xx(9 2 ) (yXi(0i) + yx 3 (0i)j = x 3 (0 2 ) + yX 3 (0 X ) 

/&- b+ \ (b- b + 

X 2 {6 2 ) f — X 2 (0x) + yX 4 (0l)j = X 4 (0 2 ) ( yX 2 (0!) + —X^) 
X 3 (6 2 ) (yX 3 (0l) + yX 4 (0!)^ = X 4 (0 2 ) ^yX 3 (0l) + y^l) 

xi(9 2 ) (x 2 {9 x ) + x 3 (9 l ) + jjx 1 (9 1 ) + ^Cx^ 



(A.6) 
(A.7) 
(A.8) 



x 4 (0 2 ) ( x 2 (6 1 ) + x 3 (0i) + z T 1 Y x 1 (e 1 ) + — ^x 4 (0i) 



d H c 
11 
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(A.9) 



x 2 (9 2 )(x 1 (9 l )+x,(9 l ) r ~ (1 ~ ' ,/_r 



+ yy^i) + jjx 3 (9 1 ) 



X 3 (6 2 ) (xi(0!) + X 4 (0!) + y yX 2 (0i) + y yX 3 (0i) 



c (i 



(A.10) 



Zi(0 2 ) ^y^i(^i) + yx 2 (0i) + yx 3 (0i) + yx 4 (0x) 
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= x 2 (e 2 ) (jrxiidt) + jx 2 (e 1 ) + jx 3 (e 1 ) + jxM^ 

X X {9 2 ) (j-X^) + y2T 2 (0i) + jX 3 (9) + yX 4 (#i)) 
= x 3 {9 2 ) ^jx 1 (9 1 ) + jX 2 (9 1 ) + LxM + y* 4 (0i)) 



(A.ll) 



(A.12) 



c~ d + \ 

f - , / + jx 3 (e 1 ) + —x 4 (e 1 )j 

(c + dr c~ d + 

Xa{0 2 ) ^-jX X {e x ) + -yX 2 {0 1 ) + yX 3 (0 X ) + yX 4 (#l) 



x 2 {6 2 ) i^jx l {9 1 ) + LxM 

dr 



\f / / / 

x 3 (6 2 ) (jx l {9 1 ) + yar a (0i) + y* 3 (0i) + y^i)) 

= X 4 (£ 2 ) (ySl(0l) + y* 2 (0l) + yX 3 (0 X ) + y3 4 (0l)) . 

Due to the relations ( |A.4j ), the functional equations (|A.11|) -( [A.14|) reduce to 
xi(6 2 ) (j-Xtfa) + jx^ij = x 2 (6 2 ) (yz 3 (0i) + y^i)) , 

xi(6 2 ) (jrx 2 {9 x ) + jx^ij = x 3 (6 2 ) (jxM + yr 4 (^ , 

x 2 {9 2 ) (y*i(0i) + jxM^ = x A {9 2 ) (jx 1 (9 1 ) + jx 3 (9x)\ , 

x 3 (9 2 ) (jx^x) + yx 2 (0i)) = x 4 (9 2 ) (yZi(0i) + yz 2 (0i)j . 

We shall solve the functional equations (|A~5| ) - (|A.10|) and ( |A.15| )- (|A.18|) . First, let 
( |A.5|) . Using ( |A.4| ), we find that ( |A.5| ) is equivalent to 



(A.13) 



( vm -..-/_ e -/n c h? x 2 (9 2 ) 
cotfr XM , q x = cot^ 2 — ^3 

From ( A.19|) , we get 



6 x 1 (9 2 ) 6 

P h 2 I h, X 2{9 2 ) 



= constant. 



ar 2 (0) _ e' h + a 1 e h tan9 

xi (9) e h — a\e~ h tan 9 ' 

where a% is a constant. Similarly, from ([A . 6|) (|A . 8|) , we get 

x 3 {9) e~ h + a 2 e h tan9 

xi(9) ~ e h - a 2 e- h ttm9' 

x 2 {9) e~ h + a 3 e h tan9 

x^{9) e h — a 3 e~ h tan 9 ' 

x 3 (9) e~ h + a 4 e h tan0 

£ 4 (0) e ft — a 4 e~ h tan ' 



(A.14) 

(A.15) 

(A.16) 

(A.17) 

(A.18) 
us consider 

(A.19) 

(A.20) 

(A.21) 
(A.22) 
(A.23) 
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where a,- 



2,3,4) are constants. Now from (|A.15 ) we have 

Jl 2 X 2 {9 2 ) _ -h 2 

= cot Q\ 



cot 9o 



6 6 xM] 



e h 2 _|_ e -ft2 



^2( 



+ e- 



(A.24) 



xi(0 2 ) ^4^1; 
Substituting ([A.20Q and ( |A.23|) into (|A.24j) , we find a constraint between the constants 

ax = -a 4 . (A. 25) 

Similarly, from (|A.16 ), we find 

a 2 = -a 3 . (A.26) 

Equations ( A.17|) and ( |A.18| ) give the constraints ( A.26 ) and (|A.25| ) respectively. Finally let us 
consider a consistency condition 



x 2 {9)x 4 {9) _ x 3 {9)x 4 {9) 
x 1 (9)x 2 (9) Xl (9)x 3 (9Y 

With the constraints (|A~25D and ( |A~26l) , ( |A~2l ) gives a relation 

e _/l + aie h tan 9 e h + a 2 e~ h tan 6> e~ h + a 2 e h tan 9 e h + aie _?t tan # 
e ft — «ie _/l tan 9 e~ h — a 2 e h tan 9 e h — a 2 e~ h tan 9 e~ h — a\e h tan 9 ' 

from which we find 

Thus there are two possibilities for Xi{9) (i — 1, • • • , 4). 

a) ai = a 2 = —a 3 = —a 4 = a 



(A.27) 



(A.28) 
(A.29) 



x 2 (9) 


e~ h + ae h tan 9 


x 3 (9) 


e~ 


h + ae h tan 9 


xi(6) 


e h — ae~ h ta.^a.9 , 


xi(0) 


e h 


— ae~ h ta.n9 


x 2 (9) 


e~ h — ae h tan# 


x 3 (9) 


e~ 


h — ae h tan 9 


x 4 (9) 


e h + ae~ h tan 9 ' 


x 4 (6) 


e h 


+ ae~ h tan 9 


= —0:4 


= a 








x 2 {9) 


e~ h + ae h t&n9 


x 3 {9) 


e~ 


h — ae h tan 9 


xi(9) 


e h — ae~ h tan#' 


xi(0) 


e h 


+ ae~ h tan 9 


x 2 (9) 


e~ h + ae ft tan# 


x 3 (9) 


e" 


h — ae h tan 9 


x 4 {6) 


e h — ae~ h tan#' 


x 4 {9) 


e h 


+ ae~ h tan 9 



(A.30) 



(A.31) 



Here a is an arbitrary constant. It is easy to confirm that both cases (|A.30|) and (|A.31|) satisfy 
the remaining functional equations ( |A.9[ ) and ( |A.10| ). To conclude, we have shown that there 
are two diagonal solutions of the RE ( |A.1| ). The solutions (|A.30|) and ( |A.31|) are equivalent to 
( |5.5| ) and (^B|), respectively. 



23 



References 

[1] E. H. Lieb and F. Y. Wu: Phys. Rev. Lett. 20 (1968) 1445. 

[2] V. E. Korepin and F. H. L. EBler: Exactly Solvable Models of Strongly Correlated Electrons 
(World Scientific, Singapore, 1994). 

[3] O. J. Heilmann and E. H. Lieb: Ann. N.Y. Acad. Sci. 172 (1971) 584. 

[4] C. N. Yang: Phys. Rev. Lett. 63 (1989) 2144. 

[5] C. N. Yang and S. C. Zhang: Mod. Phys. Lett. B 4 (1990) 759. 

[6] M. Pernici: Europhys. Lett. 12 (1990) 75. 

[7] I. Affleck: Physics, Geometry and Topology, ed. H. C. Lee (Plenum Press, New York, 1990) 
p. 1. 

[8] F. H. L. EBler, V. E. Korepin and K. Schoutens: Phys. Rev. Lett. 67 (1991) 3848. 

[9] F. H. L. EBler, V. E. Korepin and K. Schoutens: Nucl. Phys. B 372 (1992) 559. 

[10] F. H. L. EBler, V. E. Korepin and K. Schoutens: Nucl. Phys. B 384 (1992) 431. 

[11] F. H. L. EBler and V. E. Korepin: Phys. Rev. Lett. 72 (1994) 908. 

[12] F. H. L. EBler and V. E. Korepin: 
Nucl. Phys. B 426 (1994) 505. 

[13] B. M. Shastry and B. Sutherland: Phys. Rev. Lett. 65 (1990) 243. 

[14] M. J. Martins and R. M. Fye: J. Stat. Phys. 64 (1991) 271. 

[15] H. Schulz: J. Phys. C: Solid State Phys. 18 (1985) 581. 

[16] H. Asakawa and M. Suzuki: J. Phys. A: Math and Gen. 29 (1996) 225. 

[17] M. Shiroishi and M. Wadati: J. Phys. Soc. Jpn. 66 (1997) 1. 

[18] T. Deguchi and Y. Yue: preprint OCHA-PP-84 (1996). 

[19] B. S. Shastry: Phys. Rev. Lett. 56 (1986) 1529. 

[20] B. S. Shastry: Phys. Rev. Lett. 56 (1986) 2453. 

[21] B. S. Shastry: J. Stat. Phys. 50 (1988) 57. 

[22] E. Olmedilla, M. Wadati and Y. Akutsu: J. Phys. Soc. Jpn. 56 (1987) 2298. 

[23] E. Olmedilla and M. Wadati: Phys. Rev. Lett. 60 (1988) 1595. 

[24] M. Shiroishi and M. Wadati: J. Phys. Soc. Jpn. 64 (1995) 57. 



24 



[25] M. Shiroishi and M. Wadati: J. Phys. Soc. Jpn. 64 (1995) 2795. 
[26] M. Shiroishi and M. Wadati: J. Phys. Soc. Jpn. 64 (1995) 4598. 
[27] E. K. Sklyanin: J. Phys. A: Math. Gen. 21 (1988) 2375. 
[28] H. Q. Zhou: Phys. Rev. B 54 (1996) 41. 

[29] X. W. Guan, M. S. Wang and S. D. Yang: Nucl. Phys. B 485 (1997) 685. 
[30] P. P. Kulish and E. K. Sklyanin: J. Sov. Math. 19 (1982) 1596. 
[31] M. Shiroishi, H. Ujino and M. Wadati: preprint (1997). 

[32] F. Gohmann and S. Murakami: preprint ( |cond-mat / 96 10096| ) , to appear in J. Phys. A: 
Math. Gen. (1997). 

[33] L. Mezincescu and R. I. Nepomechie: Quantum Field Theory, Statistical Mechanics, Quan- 
tum Groups and Topology, eds. T. Curtright et al. (World Scientific, Singapore, 1992) p. 
200. 

[34] P. P. Kulish: Low-Dimensional Models in Statistical Physics and Quantum Field Theory, 
eds. H. Grosse et al. (Springer- Verlag, Berlin, 1996) p. 125. 



P. B. Ramos and M. J. Martins: J. Phys. A: Math. Gen. 30 (1997) L195. 
[36] R. Yue and T. Deguchi: J. Phys. A: Math. Gen. 30 (1997) 849. 



25 



